Stable magnetic monopoles in spinor polariton condensates 
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We consider half-integer topological defects in a spinor Bose-Einstein condensate of exciton- 
polaritons. We show that they accelerate under the effect of an in-plane effective magnetic field, 
like regular electric charges in the electric field. The trajectories of these analogues of magnetic 
monopoles are derived analytically and checked by numerical simulations. We show that these 
"magnetic charges" propagate at large in-plane velocities and remain stable up to a critical mag- 
netic field, thanks to the unique spin-anisotropy of polariton interactions. 

PACS numbers: 71.36.-|-c,71.35.Lk 



In the famous work of 1931 which started the 
quest for magnetic monopoles, P. Dirac considered non- 
trivial "nodal lines". These objects are characterized by 
a vanishing wavefunction on these lines and a noninte- 
grable phase around them. He found that the end of a 
nodal line is a magnetic pole with a quantized magnetic 
charge. Twenty years later, Onsagcr and Fcynman 0, Q 
described vortex lines in a 3D Bose-Einstein condensate 
(BEC) which show strong similarities with Dirac's nodal 
lines. Contrary to the nodal lines, the vortex lines termi- 
nate at the boundaries of the BEC. Another important 
difference between the nodal lines and the vortex lines 
is that the former are defined for a true single-particle 
wavefunction, whereas the latter concern the order pa- 
rameter of the condensate described with a mean-field 
Gross-Pitaevskii equation Q . Vortices are now much bet- 
ter known to scientists, because they have been nucleated 
and thoroughly analyzed experimentally (Hj , contrary to 
magnetic monopoles. 

A scalar condensate does not provide a vector field, 
which could be arranged in the required pattern, and 
therefore, an analogue of a Dirac nodal line can only be 
implemented using topological excitations in a more com- 
plex media. The similarity between the spin texture of 
certain vortices in spinor condensates and the field of a 
monopole has been rediscovered (after the initial idea of 
P. Dirac) about 10 years ago 0,01, and proposals have 
been made in order to organize the required structures 
in atomic condensates [7H9|. Moreover, the textures of 
similar type have been obtained in different media, such 
as the spin ice, and it was confirmed that they indeed be- 
have like charges, making possible what was called "mag- 
netricity" in analogy with "electricity" [lol |. 

In this Letter, we show the striking analogy between 
magnetic monopoles and half-integer topological defects 
in spinor polariton condensates. Polaritons are the quasi- 
particles formed of photons and excitons strongly coupled 
in microcavities [ll| . They have recently become a model 
system for studying condensation |l2| and various asso- 
ciated effects such as superfluidity [13| , formation Il4l 
and evolution [lH of vortices, oblique solitons [3, Il7l |. 



and others. Their popularity is backed up by their fas- 
cinating properties, mixing those of light (small effective 
mass) and matter (self-interactions and thermalization 
with phonons). Their spin structure is especially interest- 
ing: polaritons are bosons with only two allowed spin pro- 
jections on the growth axis: ±1. The most important fea- 
ture of polaritons for the present work is the strong spin- 
anisotropy of their interactions [isj . Indeed, the exciton- 
cxciton interaction constant in the triplet configuration 
ai is repulsive and approximately ten times larger than 
in the singlet configuration a2 [l9[ . Consequently, the in- 
teraction energy of the polariton condensate is minimized 
when the latter is linearly polarized p^: n+ = ri- = n/2 
(n is the total density). The spin-anisotropy of interac- 
tions makes the condensate stable against small in-plane 
effective magnetic fields which are acting on the polariton 
pseudospin |18j, as we will see below. 

The two spin components make the existence of half- 
vortices 21 1 and half-solitons(20j in polariton condensates 
quite natural. Indeed, the first experimental observa- 
tion of a half- vortex has been carried out for a polariton 
condensate |lj|. We demonstrate that a constant effec- 
tive magnetic field (whose origin will be described be- 
low) accelerates such half-integer defects, which there- 
fore behave as magnetic charges capable of propagating 
at high velocities. We demonstrate that these defects re- 
main stable if the field does not exceed a critical value. 
The ID half-soliton case is studied first, followed by a 
more complicated configuration of a half-vortex in a 2D 
condensate. 

ID system Let us consider a two component conden- 
sate described by the spinor Gross-Pitaevskii equations 
(CPE): 



in 



dt 



2m 



-Al/^cr + ailV'o-I^V'cr - TT'^-o-, (1) 



where ct = ± and m* is an effective mass associated with 
a parabolic dispersion approximation for polaritons, valid 
at relatively small wave vectors. The term with VL de- 
scribes the linear polarization splitting or in-plane effec- 
tive magnetic field, mixing the circular components. The 
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FIG. 1. (Color online) a) Wavefunction components for a dark 
half-soliton together with the corresponding in-plane pseu- 
dospin pattern; b) Trajectory of a half-soliton accelerated in 
a constant in-plane effective magnetic held calculated analyt- 
ically (red solid line) and numerically (black points). 



spin anisotropy of the polariton-polariton interaction is 
fully taken into account by using the conditions ai ^ 02 
and neglecting a2 hereafter. For the follovifing discussion, 
it is convenient to use the pseudospin representation : 

Sy,s,) with = 3?(V'+^l), Sy = 



(5, 



s = , 

Sz = ('^+ — "■-) /2. Indeed, the pseudospin evolution 
can be mathematically described by the interaction of 
this pseudospin with an effective magnetic field, similarly 
to that of an ordinary magnetic momentum subject to a 
real magnetic field 18|. The in-planc effective magnetic 
field in microcavities can describe for example the energy 
splitting between the two orthogonal linear polarizations 
linked with crystallographic axes and it can be controlled 
by an applied electric field There is usually an ad- 
ditional LT splitting in ID wire cavities [24 1. 

For a homogeneous solution -0^ = n/2, ipzfi — 
y^n/2e^'^° , the relative phase 4>o between the two com- 
ponents gives the direction of the linear polarization, or 
the in-plane component of the pseudospin S|| = {Sx, Sy). 
It is well known that the GPE admits soliton solu- 
tions, the latter can lie in only one of the two compo- 
nents of a spinor condensate giving a half-soliton (HS): 
/7T,/2tanh {x/£^\/2) where ^ 



V'+(a;) = ■\/n/2tanh (a;/^v2) where ^ ~ h/^ainm* is 
the healing length of the condensate, while the other 
component can remain constant ip- = ■\/n/2e*'^". The 
in-planc pseudospin in the particular case where 0o = 
is given by Sx{x) ~ ntanh (a;/^-\/2), Sy{x) — 0. Far 
from the soliton's core of the width ^, the pseudospin 
has therefore exactly the pattern of the field created by 
a single magnetic charge in ID, shown in Fig[lja). 

Let us analyze in the limit L ^ ^, the effect of a mag- 
netic field on the HS in order to verify that it behaves 
indeed as a magnetic charge. To find the force acting 
on the HS we have to find the gradient of the magnetic 
energy of the system as a function of the displacement of 
the soliton's core. The magnetic energy of a condensate 
containing a HS in an external in-plane magnetic field is 
/ (r2/2 • S)dx, with S(x — .to) = nsign(a; — xo)ex in the 
limit we consider, giving E{x) ~ — fi^^nx and F = Q.xn. 
For a grey soliton propagating with a speed v, the phase 



shift induced by the soliton is A5' = 2arccos(w/c) < tt 
and the pseudospin projection Sx is reduced, which can 
be expressed as a renormalization of the magnetic charge. 
In the limit L ^ the correction to the charge is found 
analytically as q = go(l ~ u^/c^). The total correction 
for the mass of the soliton and its charge gives the equa- 
tion of motion a = qaVt^n / ma{l — w^/c^)^/^, the same 
as in relativistic physics, integrating which one obtains 
v{t) = cta,Vi\i{qoQ.x'nt / c) (zero initial velocity). The re- 
sulting trajectory is shown in Figlljb), perfectly fitting 
the results of numerical simulations (see below) . 

Under the effect of the effective magnetic field the HS 
can indeed propagate similarly to a real charge, its stabil- 
ity against the pseudospin rotation being maintained by 
the interaction energy. However, if the in-plane magnetic 
field exceeds a certain value, the HS can change sign: 
convert into the other polarization component. The esti- 
mation of the critical field J7c can be obtained analyzing 
the dynamical equation for the pseudospin in the center 
of the HS: dtS = Jl x S, with = Qx/fi^x + aiS'^e^. 
The analytical solution of this equation is expressed in 
terms of elliptical functions demonstrating regular oscil- 
lations, and the asymptotic behavior of the S^ pseudospin 
component at critical field fic — ain^/2 is exponential. 
Below flc the pseudospin keeps the initial sign of its pro- 
jection over z: the HS is stable. In atomic condensates 
the spin anisotropy is usually almost absent (ai ~ 0:2) 
and thus the HS would be destroyed. The same result 
applies to the 2D case for half-vortices as we will see. 
This dynamical nonlinear effect is different from the spin- 
Meissner effect obtained for the ground state of polari- 
tonic systems [2^. 

To check the analytical predictions, we have performed 
two types of numerical simulations using the spinor GPE. 
First, we have imposed a HS as an initial condition and 
simulated its propagation under the effect of a magnetic 
field. The resulting circular polarization degree is shown 
in Figl2][a), the blue region corresponding to the soliton. 
The time dependent position of the soliton has been ex- 
tracted from this simulation and successfully compared 
with analytical predictions in Fig[ljb) . 

Second, we have performed a more realistic simulation 
of a wire-shaped cavity with finite lifetime under quasi- 
resonant, homogeneous, and linearly polarized pumping 
by an external laser of amplitude and frequency ui: 

(2) 

In this case the HS is created after the steady state 
regime is obtained, by a short potential pulse acting on a 
single polarization component cr"*". Applying this poten- 
tial on the wire's edge allows creating a single HS propa- 
gating at speed v and accelerating as its depth decreases. 
Thus, the finite lifetime acts as an effective force [1], ac- 
celerating a relativistic particle (m ~ mg / \/l — j (?) 
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FIG. 2. (Color online). Circular polarization degree as a func- 
tion of coordinate and time calculated using the spinor GPE, 
including the effect of a constant in-plane effective magnetic 
field. The HS trajectory is visible as the deep blue minimum, 
a) initial conditions, no lifetime; b) pumping and lifetime, the 
HS is created by a pulsed potential. Black dashed lines are 
the guides for the eyes showing trajectories for the splitting 
values indicated on the figure (in meV). 



up to the limiting speed, proportional to the speed of 
sound c = ^ ain/2m. The effect of the magnetic field 
can be seen as an extra force F{t) ~ q(v{t))Qxn, which 
can decrease or increase the acceleration [Figl2]^b)]. We 
see that even in the driven-dissipative case the magnetic 
monopoles are suflaciently stable to propagate under the 
effect of the magnetic field on long distances. 

An important difference between the Dirac's 
monopoles and the half-integer defects in spinor 
condensates lies in the interactions between them. 
While a real magnetic charge creates a real magnetic 
field which certainly influences other magnetic charges, 
half-integer defects possess only a (pseudo)spin field, 
which has a weak influence on other half-integer defects. 
They behave indeed as magnetic charges in an external 
magnetic fleld, but see each other rather via interaction 
and kinetic energy. 

2D system In the 2D case the object possessing the 
spin texture close to that of a magnetic monopole is the 
elementary half-vortex (HV) carrying winding number 
±1 in the component where it appears. However, there 
is an important difference between the fleld of an electric 
charge, which in 2D decreases as 1/r, and the pseudospin 
"fleld" of the HV, which is approximatively constant at 
large distances. In spite of that, it is still possible to find 
the force acting on the HV from the magnetic field as the 
gradient of the energy of the system as the function of 
the displacement of the monopole. 

A crucial peculiarity arises from the relative phase 
(l)Q of the two wavefunction components: ■;/'+ = 



/Ti-e 



Ki-'t>) 



is the polar angle) 



which gives an extra degree of freedom for the pseudospin 
orientation. Several pseudospin textures are shown in 
FiglSl^a) together with the resulting sketched trajectory 
of the HV in the constant in-plane magnetic field 
(black arrow). These trajectories have been calculated 
by solving the 2D spinor GP equations with infinite life- 
time. We have considered a cylindrical trap of radius 



50 ^m with impenetrable boundaries and we have found 
the associated HV "ground state" at fixed chemical po- 
tential /i = 1 meV with a given set (/+, 0o). Then, for 
each of the 8 HV configurations shown on the Fig. 3(a) 
we have applied a constant magnetic field along the x- 
direction. The resulting calculated trajectories are shown 
in black on the Fig. 3(b). Experimentally, a half-vortex 
can be created in a cylindrical trap (mesa) with a Gauss- 
Laguerre laser beam. 
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FIG. 3. (Color online), a) Pseudospin textures for HVs as a 
function of topological charges and relative phase cj>o (labeled 
as (/+,?_, (^o)) and the resulting propagation directions in a 
constant magnetic field pointing along the x direction, b) 
The (j^ density at the end of the simulation with the Gross- 
Pitaevskii equation showing the trajectory (red line) of a par- 
ticular HV (white square) together with the other trajectories 
(black lines) corresponding to the cases described on a). 

We see that the propagation direction is determined 
both by the topological charges (/+,/_) of the HV and by 
the continuous variable c^o- We note that the (Z_|_, Z_, (/ig) 
and the (— — L, (/)o) HVs are symmetric with respect 
to the magnetic field and follow the same trajectory. 
Thus, the proper description for the force acting on the 
HV from the magnetic field is given by a charge tensor: 
Fi = QijHj, where qij are proportional to components of 
the rotation matrix for angle (tt— (^o) (for the conventions 
used above). 

In order to understand better the forces acting on the 
HV, let us consider the (1,0,0) HV, whose trajectory is 
shown in red and whose pro pag ation is demonstrated in 
the supplementary material |27| . Two effects induced by 
magnetic field, are observed in the simulations. 

First, because of the pseudospin rotation around the 
magnetic field far from the vortex core, a density gradi- 
ent appears at short times along the field's axis {x in our 
case) as shown on the [Fig|4|a)]. The density determines 
the vortex energy, whose gradient creates an additional 
force acting on the HV. Fortunately, this force is strong 
and well defined only in the initial moments after the 
application of the field, because later on the different 
eigenfunctions of the vessel oscillate with different fre- 
quencies and the resulting " storm in a teacup" gives zero 
net force for the vortex, which therefore propagates with 
a constant speed acquired in the initial moments. 



4 



Second, the constant effect of the magnetic field on the 
pseudospin texture creates a constant force accelerating 
the vortex in the direction given by eq. (3). The tra- 
jectory is parabolic, with the acceleration proportional 
to the strength of the magnetic field, as shown on the 
[FigHUb)]. It is precisely this effect, and not the pre- 
vious one, which corresponds to the expected behavior 
of a magnetic monopole. The HV is remarkably sta- 
ble against density fluctuations and we are able to track 
its trajectory for 80 ps. Relativistic effects are not ob- 
served in this case, because the speed remains rela- 
tively small. 

We note that the effect of a static polarization splitting 
on a HV has been analyzed in Ref. |26| and a polarization 
string has been predicted to be attached to the HV in the 
lowest energy state. This is indeed what we observe in 
the early moments of the simulation, when the magnetic 
field is turned on [27| , but the string tends to unbind from 
the HV because our initial condition corresponds to the 
absence of the field and this initial condition is not the 
eigenfunction of the complete Hamiltonian, otherwise it 
would not evolve with time. 
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FIG. 4. (Color online), a) The density distributions of the two 
components: (red) and cr~ (blue), creating the gradient of 
the interaction energy accelerating the HV in the y direction, 
b) The fitted a;-coordinate of the HV as a function of time for 
several values of the magnetic field, demonstrating constant 
acceleration of the half-vortex, increasing with the field. 

The HV can become unstable against magnetic field, 
with the stability criterion given in our case by the same 
estimate as in ID. If the magnetic field is stronger, the 
vortex disappears. Without spin-anisotropy of interac- 
tions in the homogeneous magnetic field we consider, the 
pseudospin along the y axis would make a turn around 
the magnetic field and become completely Sz on one 
side and —Sz on the other side, leading to the destruction 
of the HV. Roughly speaking, instead of having a single 
point with zero density of ct+, we would have a semi- 
infinite string. This is the signature of the fact that in 
the case of a spinor condensate, the velocity circulation 
is not conserved topologically, because one cannot define 
unique superfluid velocity |28|. Therefore, the HVs can 
be destroyed (unwound) by applying a certain magnetic 
field, as it has been demonstrated experimentally [29[. In 
the case of polaritons, the particular spin-dependence of 
interactions prevents this from happening for sufficiently 



low magnetic fields. 

Finally, we would like to note that the relative phase 
between the condensate components can be modified by 
applying a weak external magnetic field in the z direction, 
leading to a Zeeman splitting for the (t"*" and a~ compo- 
nents. The relative phase would increase linearly with 
time and therefore the direction of the in-plane magnetic 
force would change as well. 

Conclusions Half-integer topological defects in spinor 
polariton condensates behave as magnetic charges with 
respect to an in-plane effective magnetic field. Their 
stability is provided by the spin-anisotropy of polariton- 
polariton interactions. We have shown that the dynam- 
ics of half-solitons in ID is well described by relativistic 
equations. 
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